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SUMMARY 

Reliable  experimental  data  on  local  heat-transfer  coefficients  for 
supersonic  flow  of  air  in  a round  tube  are  reanalyzed  in  detail  w’ith  the  aid 
of  an  approximate  two-dimensional  flow  model.  The  results  are  compared 
with  similar  results  based  on  a one -dimensional  flow  model  and  with  the 
theoretical  predictions  for  supersonic  flow  over  a flat  plate  and  for  flow  in 
the  entrance  region  of  a tube  when  a laminar  boundary  layer  is  present. 

The  two-dimensional  flow  model  yields  a better  understanding  of  the 
phenomena  which  occur  for  diabatic  supersonic  flow  of  air  in  a round  tube 
than  that  obtained  with  the  aid  of  the  one-dimensional  flow  model.  The  two- 
dimensional  flow  model  shows  that  the  core  Mach  number  is  nearly  constant 
along  the  length  of  test  section  for  a range  of  values  of  the  inlet  diameter 
Reynolds  number.  For  a laminar  boundary  layer,  the  values  of  the  local 
Stanton  number-  agree  within  a few  per  cent  with  the  theoretical  values  for 
p'ate  flow  at  the  largest  values  of  the  inlet  diameter  Reynolds  number. 
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NOMENCLATURE 


a - inside  radius  of  pipe 

A - cross-sectional  area,  7rD“/4 

A!  - heat -transfer  area,  tUAL 

Cp  - specific  heat  constant  pressure 

cw  - discharge  coefficient  of  nozzle 
D - inside  diameter  of  pipe 

g - acceleration  given  to  unit  mass  by  unit  force 
G - flow  per  unit  area,  w/A 
h - coefficient  of  heei,  transfer,  4/ A1  (tw-taw) 
k - ratio  of  specific  heats 

L - distance  from  end  of  curved  contour  of  nozzle 
M - Mach  number,  V I V gkRT 
n - summation  index,  Eq.  (15) 

Nul  - length  Nusselt  number,  iiL/A 

p - static  pressure 

q - net  rate  of  heat  transfer 

R - perfect  -gas  constant 

Rejy  - diameter  Reynolds  number,  DG/pg 

ReL  - length  Reynolds  number,  LG/pg 

St  - Stanton  number,  h/c  G 

t - temperature,  deg  F 

T - temperature,  deg  F abs 

V - velocity 

w - mass  rate  cf  flow 

p - density 

P - viscosity 

A - thermal  conductivity 

5 - thickness  of  boundary  layer 

Superscript  * refers  to  throat  of  supersonic  nozzle  where  M = 1. 


Subscripts: 

aw  - adiabatic  wall  conditions 
b - boundary  layer 

c - isentiopic  core 

j - station  numbers 

0 - hypothetical  entrance  plane  of  the  tube,  where  the  boundary  layer  is 

of  zero  thickness 

ob  - stagnation  conditions  in  boundary  layer 

oc  - stagnation  conditions  in  core 

01  - upstream  stagnation  conditions 

r - atmospheric  conditions 

s - isentropic  conditions 

w - waR  conditions 

00  - free  stream  conditions  for  flat-plate  flow 


INTRODUCTION 


An  inexpensive  steady-state  supersonic  wind  tunnel  has  been  designea 
Tor  accurate  measurements  of  'local  heat -transfer  coefficients.  It  consists  of 
a convergent-divergent  nozzle  which  delivers  dry  air  at  supersonic  speeds  to 
a test  section,  one-half  inch  in  diameter  and  30  to  50  diameters  in  length.  The 
various  test  apparatus  and  the  measurements  obtained  therefrom  have  been 
described  previously  in  detail  (1).*  Sufficient  measurements  have  been  made 
to  demonstrate  their  reliability  when  a laminar  boundary  layer  exists  in  the 
test  section  for  adiabatic  and  diabatic  supersonic  flow.  The  heat -transfer  data 
were  first  analyzed  (1)  with  the  aid  of  a simple  one-dimensional  flow  model, 
abbreviated  hereafter  as  1 -DFM,  but  it  was  shown  that  this  model  is  not  fully 
adequate  to  explain,  interpret  and  correlate  these  measurements. 

The  value  of  these  measurements  of  local  heat-transfer  coefficients  would 
be  considerably  enhanced  if  they  could  be  compared  with  similar  measurements 
for  other  types  of  supersonic  flew  with  a laminar  boundary  layer,  such  as  plate 
flow,  cone  flow,  etc.  One  reason  for  subjecting  the  data  given  in  (1'  to  further 
analysis  and  computation  is  to  lay  the  foundation  for  such  comparisons.  A 
second  reason  is  that  accurately  measured  values  of  local  heat-transfer  co- 
efficients for  supersonic  flow  with  a laminar  boundary  layer  are  practically 
nonexistent  so  that  further  analysis  of  tube -flow  data  might  produce  informa- 
tion useful  to  the  de signet*  oi  devices  moving  at  supersonic  speeds.  A third 
reason  for  using  a more  complicated  flow  model  to  reanalyze  the  same  data  is 
to  provide  the  exnerimenter  with  a better  phenomenological  explanation  of  the 
processes  which  occur  in  supersonic  now  in  ihe  entrance  region  of  a tube.  The 
two-dimensional  flow  model,  abbreviated  hereafter  as  2-DFM,  will  be  used  in 
this  paper  to  analyze  and  interpret  some  selected  values  of  the  data  given  in(l). 

The  analysis  of  supersonic  flow  in  the  inlet  region  of  a tube  can  be  ac- 
complished, on  the  one  hand,  in  an  extremely  simple  fashion,  such  as  with  the 
aid  of  ihe  1-DFM,  and  on  the  other  hand,  in  an  extremely  difficult  fashion, 
such  as  by  investigation  of  the  partial  differential  equations  of  energy,  momen- 
tum, and  continuity.  Although  both  of  these  methods  ol  analysis  have  been  used 
in  the  present  research  program,  an  intermediate  type  of  analysis  is  used  here 
to  obtain  results  close  to  those  based  on  the  most  exact  analysis.  This  2-DFM 
should  be  considered  to  represent  closely  the  actual  supersonic  flow  phenomena 
in  the  tube. 

The  2-DFM  is  used  here  only  for  supersonic  flew  in  the  entrance  region 
of  a tube  with  a laminar  boundary  layer  originating  at  tube  inlet.  Experimentally 
it  has  been  found  possible  to  maintain  a laminar  boundary  layer  over  the  entire 
length  of  test  section  by  careful  control  of  the  inlet  conditions.  Independent 
measurements  of  the  velocity  and  temperature  profiles  of  the  tube  flow  for 
these  inlet  conditions  have  confirmed  the  existence  of  this  laminar  boundary 
boundary  layer.  In  the  previous  paper  (1)  ceriain  values  of  the  tube  inlet  con- 
ditions were  selected  to  present  the  data  for  a laminar  boundary  layer.  In 

♦Numbers  in  parentheses  refer  to  the  Bibliography  at  the  end  of  the  paper. 


2 


the  present  paper,  the  same  values  of  these  inlet  conditions  are  selected  to 
present  the  results  of  computation  based  on  the  2-DFM.  Comparisons  of  the 
results  based  on  the  1-DFM  and  on  the  2-DFM,  for  the  same  original  date  will 
be  given  here.  The  results  based  on  the  2-DFM  will  also  be  compared  with 
theoretical  solutions  for  plate  flow  and  for  tube  flow. 


ASSUMPTIONS  FOR  TWO-DIMENSIONAL  FLOW  MODEL 
WITH  HEAT  TRANSFER 

The  2-DFM  arbitrarily  divides  the  supersonic  flow  in  the  entrance  region 
of  a tube  into  two  parts,  a laminar  boundary  layer  and  a central  core  of  fluid, 
as  shown  in  Fig.  1.  The  laminar  boundary  layer  originates  near  the  tube  inlet 
and  grows  in  thickness  by  transference  of  mass  from  the  central  core.  This 
2-DFM  is  consistent  with  Prandtl1  s concept  of  the  boundary  layer  since  it  is 
assumed  that  both  the  viscous  and  thermal  effects  are  concentrated  in  the 
boundary  layer,  and  that  both  are  negligible  in  the  core  region,  i.e.,  the  fluid 
remaining  in  the  central  core  undergoes  an  i3entropic  process.  The  basic 
soundness  of  this  2-DFM  has  been  conclusively  demonstrated  by  means  of 
accurately  measured  velocity  profiles  at  various  positions  along  the  test  section 
of  Fig.  1,  for  both  adiabatic  and  diabaiic  flow  (2,  2),  and  also  by  means  of  un- 
published recentl>  measured  temperature  profiles. 

The  present  paper  deals  with  a simplified  version  of  the  actual  2-DFM  in 
that  the  actual  velocity  and  temperature  profiles  are  represented  by  simple 
approximations  for  both  the  laminar  boundary  layer  and  the  core.  The  boundary- 
layer  thickness  is  assumed  to  be  the  same  for  both  the  velocity  and  thermal 
houndary  layers.  This  thickness  is  arbitrarily  taken  to  correspond  to  that 
thickness  at  which  the  local  velocity  in  the  boundary  layer  equals  99  per  cent 
of  the  mean  velocity  in  the  core.  These  simplifications  are  shown  in  Fig.  1. 

They  were  introduced  mainly  to  reduce  the  computational  time  since  i+  was 
found,  after  considerable  study,  that  the  use  of  these  simplified  velocity  and 
temperature  profiles  yielded  results  nearly  identical  to  those  obtained  using 
a non-linear  velocity  profile  but  requiring  considerably  greater  computational 
time. 


The  following  assumptions  are  made  in  the  analysis  of  the  2-DFM: 

1.  A laminar  boundary  layer,  in  which  both  the  viscous  and  thermal 
effects  are  predominant,  exists  for  supersonic  flow  in  the  entrance  region. 

2.  The  air  which  remains  in  the  central  core  undergoes  an  isentropic 
change  of  state  from  an  upstream  stagnation  state  to  the  static  pressure 
measured  at  a given  section. 

3.  The  properties  in  the  central  core  of  fluid  are  assumed  to  be  uniform 
at  a given  section. 
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4.  The  static  pressure  Is  assumed  to  be  uniform  across  each  section 
of  the  tube. 

5.  Air  is  a perfect  fjas  with  a constant  value  of  the  ratio  of  specific 
heats  (k  = 1.40)  over  the  range  of  temperature  under  consideration. 

6.  The  laminar  boundary  layer  has  a linear  distribution  of  velocity 
which  is  further  approximated  by  an  average  velocity  in  the  boundary  layer 
equal  to  one-half  the  core  velocity,  as  shown  in  Fig.  1. 

7.  The  stagnation  temperature  of  the  central  core  is  assumed  to  be 
constant  at  all  sections,  equal  to  Toc,  whereas  the  stagnation  temperature  of 
the  boundary  layer  is  assumed  to  be  uniform  at  a given  section  but  variable 
from  section  to  section. 

8.  The  flow  is  adiabatic  up  to  the  tube  inlet  where  the  heat  transfer 
originates. 


ANA. LYSIS  OF  SIMPLIFIED  TWO-DIMENSIONAL  FLOW  MODEL: 

WITH  HEAT  TRANSFER 


The  following  relations  hold  at  each  section  of  the  flow  in  the  tube  shown 
in  Fig.  1 . 


Continuity: 

w 

a wb  + wc 

(i) 

C ontinuity: 

wc 

= PrVcAc 

(2) 

Continuity: 

wb 

r PbvbAb  = Pbvb(A-Ac> 

(3) 

Geometry: 

Ac 

= * [ (D/2)  - 6^  2 

(4) 

Geometry: 

A 

= Ab  + Ac 

(5) 

Equation  of  state; 

P 

= pcRTc 

(6) 

Equation  of  state: 

P 

" PkRTb 

(7) 

Definition: 

M2 

c 

= V2/gkRT 

ft  c 

(8) 

Assumption: 

Vb 

= Vc/2 

(9) 

The  following  relations  hold  between  the  upstream  stagnation  state  and 
a state  in  the  core  at  any  tube  section: 


Isentropic: 

Isentropic: 

Isentropic: 

Isentropic: 


Pol1?  " <*oi'Vk 
Toi/Tc=  (poi/pl^-D/k 

T . = T 

Ol  oc 

c T = c T + V2/ 2g 

P oc  pc  c 


(10) 

(ID 

(12) 


(13) 


Definition: 


c p 1 ob  _ cp^b  + ^'b^S 


(14) 


Application  of  the  energy  equation  for  steady  flow  to  a control  volume 
which  encloses  the  fluid  between  the  upstream  stagnation  region  and  any  tube 
section  j yields 


(1  / 2 ) 


n = j - 1 

+ n 


n = 1 


‘n 


g 

w . (c  T . + V“/2g)  + w 
c;-  P cj  c 


bj 


(c  T,  . + V / 2g)-wc  T ..  (15) 
P bj  b p oi 


The  summation  in  Equation  (15)  starts  at  the  first  station  at- which  heat 
transfer  data  were  measured.  This  first  station  for  test  combinations  C and 
D is  described  in  reference  (1). 

The  additional  relations  used  in  this  analysis  are  as  follows:  The  condi- 
tion for  choked  flow  of  air  through  the  supersonic  nozzle  yields 


G* 

s 

= 

(w/A*)  - 0.5318  p ./  ./ T~^ 

* oi  V oi 

(16) 

Definition: 

c 

-- 

(w/A*) /(w/A*) 

(17) 

w 

s 

Assumption: 

T IT  . 

= 

0.940 

(18) 

aw  oi 

Definition: 

h 

= 

q/A'  (t  -t  ) 

(19) 

w aw 

Definition: 

NuDc 

- 

hD/X 

c 

(20) 

Definition: 

Nu, 

Lc 

- 

hL/X 

c 

(21) 

Definition: 

st. 

= 

h/c  p V 

(22) 

c 

pc  c 

Definition: 

K*Dc 

- 

p V D / ju  g 
c c c c 

(23) 

Definition: 

R6Lc 

- 

p V L/u  g 
c c c 

(24) 

Method  of  Computation 

A sample  calculation  is  given  in  the  Appendix. 


5 


RESULTS 


General 

The  original  measurements  of  heat-transfer  coeffir  ients  for  supersunic 
flow  obtained  with  test  combinations  C and  D,  comprising  four  runs  with  the 
former  and  thirteen  runs  with  the  latter,  are  given  in  detail  in  (1).  These  data 
have  been  selected  for  those  stagnation  conditions  where  a laminar  boundary 
layer  is  present  over  most  of  the  test  section.  Hence,  the  method  of  computa- 
tion outlined  above  for  the  simplified  2-DFM  should  be  applicable  for  these 
seventeen  runs.  In  view  of  the  large  amount  of  calculation  involved,  however, 
only  seven  runs  have  been  studied  in  detail  by  means  of  the  simplified  2- ITEM, 
namely  six  and  one  from  test  combinations  C and  D,  respectively. 

The  results  for  six  runs  made  with  test  combination  D are  given  in  Figs. 

2 to  7,  inclusive.  Six  charts  were  found  necessary  to  present  these  results  in 
sufficient  detail  because  they  are  sensitive  to  small  changes  in  the  value  of  the 
inlet  diameter  Reynolds  number.  The  results  for  one  run  made  with  test  com- 
bination C are  given  in  Fig.  8.  One  run  was  selected  here  mainly  to  reduce 
the  total  computational  time  and  to  show  the  effect  of  a shorter  test  section 
on  the  data.  On  each  of  these  seven  charts  the  calculated  quantities  based  on 
one  flow  model  are  plotted  without  reference  to  the  quantities  calculated  from 
the  other  flow  model.  Thus,  in  the  curve  of  Mach  number  versus  length 
Reynolds  number,  the  values  of  Mach  number  and  Reynolds  number  are  taken 
from  either  the  l-DFM  or  from  the  2-DFM. 

Each  of  the  seven  charts  presents  the  Mach  number,  bound  ary -layer 
thickness  ratio,  local  length  Nusselt  number,  and  local  Stanton  number.  For 
the  2-DFM,  the  Mach  number,  Nusselt  number,  and  Stanton  number  are  based 
on  properties  of  the  central  core  of  fluid.  The  local  Nusselt  and  Sianton  numbers 
are,  in  reality,  average  values  over  the  short  length  of  flow  of  two  or  four  pipe 
diameters.  In  the  discussion  of  these  seven  charts,  the  value  of  the  inlet 
diameter  Reynolds  number  chosen  for  identification  corresponds  to  that 
computed  by  means  of  the  2-DFM 

Each  of  the  seven  charts  contains  a comparison  of  the  results  based  on 
the  2-DFM  with  the  theoretical  predictions  of  supersonic  flow  with  a laminar 
boundary  layer  on  a flat  plate  and  with  the  theoretical  predictions  for  tube 
flow  for  constant  values  of  the  viscosi'ty  and  the  thermal  conductivity  . The 
former  predictions  are  taken  from  van  Driest  (4)  and  the  latter  from  Toong  (5). 
This  type  of  comparison  of  results  based  on  the  2-DFM  with  predictions  for 
plate  flow  and  for  tube  flow  should  be  regarded  as  a temporary  expedient  until 
exact  solutions  for  tube  flow  and  accurate  experimental  data  for  plate  flow 
with  adverse  pressure  gradients  are  available. 


Results  for  Test  Combination  D-Mach  Number 

Fig.  2 presents  the  results  based  on  the  simplified  2-DFM  for  the  lov/est 
value  of  the  inlet  diameter  Reynolds  number,  namely,  33,000.  The  Mach  number 
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based  on  the  2-DFM  decreased  only  from  2.7  to  2.).  along  the  length  ot  the  test 
section  whereas  that  for  the  1-DFM  decreases  from  2.1  to  about  1.2.  This 
same  type  of  behaviour  of  the  Mach  numbers  for  the  two  flow  models  is  found 
also  in  Figs.  3,  4,  5,  6,  and  7.  In  Fig.  7,  the  Mach  number  based  on  the  2-DFM 
decreased  only  from  2.65  to  2.50  over  50  diameters  of  length  of  flow.  This 
constancy  of  the  Mach  number  based  on  the  2-DFM  supports  the  contention 
that  this  type  of  test  section  is,  in  reality,  a useful  and  inexpensive  type  of 
steady-state  supersonic  wind  tunnel.  The  2-DFM  shows  that  the  core  Mach 
number  does  not  vary  as  much  as  that  indicated  by  the  1-DFM  and  that  ;he 
value  of  the  core  Mach  number  is  significantly  greater  than  the  corresponding 
value  based  on  the  1-DFM.  Independent  measurements  of  velocity  and  tempera- 
ture profiles  at  several  stations  near  the  tube  exit  lead  to  calculated  values  of 
the  core  Mach  number  which  agree,  within  one  per  cent,  with  the  values  of  the 
core  Mach  number  based  on  the  2-DFM. 


Boundary -Layer  Thickness 

The  boundary-layer  thickness  at  any  section  in  the  tube  is  defined  by 
Equation  (4)  for  the  2-DFM.  It  should  be  noted  that  this  thickness  is  taken 
arbitrarily  as  zero  in  the  1-DFM.  For  the  smaller  values  of  the  inlet  diameter 
Reynolds  number,  namely  from  33,000  to  72,000  in  Figs.  2,  3,  and  4,  this  ratio 
increased  from  about  20  per  cent  at  tube  inlet  to  about  40  to  50  per  cent  at 
tube  exit.  For  the  larger  values  of  Reynolds  number,  namely  from  88.000  to 
115,000  in  Figs.  5,  6,  and  7,  this  ratio  increases  from  about  10  per  cent  at  tube 
inlet  to  about  25  per  cent  at  tube  exit.  Such  behaviour  of  the  thickness  ratio 
could  be  interpreted  in  two  ways.  Either  the  laminar  boundary  layer  does  not 
fill  the  tube  at  exit  and  a fully  developed  flow  pattern  is  not  attained  for  this 
entrance  type  of  flow,  or  for  the  flow  length  of  50  diameters  the  supersonic 
flow  is,  in  reality,  one  corresponding  to  a fully  developed  flow  but  with  a velocity 
profile  which  is  quite  different  than  that  for  fully  developed  incompressible  flow 
in  a round  tube.  Additional  experimental  data  on  velocity  profiles  for  several 
stations  near  the  tube  exit  are  needed  before  these  views  can  be  resolved 

The  thickness  ratio  for  tube  flow  based  on  the  2-DFM  can  be  compared 
with  a similar  ratio  derived  for  the  thickness  of  a laminar  boundary  layer  of 
a compressible  fluid  flowing  on  a flat  plate  for  zero  pressure  gradient.  The 
theoretical  results  of  Howarth  (6)  lead  to  the  following  equation  for  the  thickness 
ratio, 

(6/a)  = 10.4(L/D)(1  + 0.0795  M2)/(Re„  )1'2.  (25) 

c Lc 

Here  the  thickr  .s  corresponds  to  that  value  at  which  the  local  velocity  in 
the  boundary  layer  on  the  flat  olate  is  99  per  cent  of  the  free  stream  value. 

The  values  of  the  thickness  ratio  based  on  the  2-DFM  arc  slightly  larger  than 
those  based  on  Equation  (2  5)  for  values  of  the  diameter  Reynolds  number  from 
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33,000  to  72,000,  but  both  sets  of  values  are  in  excellent  agreement  over  the 
entire  length  of  the  test  section  for  Reynolds  numbers  from  88,000  to  115,000. 
Tnis  excellent  agreement  should  be  considered  as  somewhat  coincidental  for 
the  following  reasons: 

1.  The  zero  of  the  length  measurement  for  the  tube  flow  should  be  the 
origin  of  ihe  boundary  layer  where  its  thickness  is  zero,  but  the  aetual  zero 
has  been  arbitrarily  taken  to  be  at  the  end  of  the  curved  contour  of  the  nozzle. 
The  effects  introduced  by  this  error  should  disappear  by  the  middle  of  the 
test  section. 

2.  The  Howarth  formula  is  restricted  to  zero  pressure  gradient  where- 
as a finite  adverse  pressure  gradient  exists  in  the  tube  flow.  The  smallest 
pressure  gradient  is  found,  however,  at  the  largest  value  of  the  inlet  diameter 
Reynolds  number,  and  for  this  case,  the  agreement  between  thickness  ratio 
based  on  the  2-DFM  and  on  the  Howarth  formula  is  the  best. 

3.  At  the  lowest  values  of  the  diameter  Reynolds  number,  namely  at 
33,000  to  72,000,  some  evidence  exists  that  separation  of  the  laminar  boundary 
layer  probably  occurs  in  the  presence  of  the  adverse  pressure  gradient.  Such 
separation  would  probably  result  in  a thickening  of  the  laminar  boundary  layer 
and  cause  its  calculated  thickness  based  on  the  2-DFM  to  be  greater  than  that 
indicated  by  the  Howarth  formula.  Such  a thickening  calculated  on  the  basis 

of  the  2-DFM,  ean  be  seen  in  Figs.  2,  3,  and  4. 


Stanton  Number 

The  local  heat-transfer  coefficients  based  on  the  2-DFM  are  shown  in 
Figs.  2 to  7 for  test  combination  D in  the  form  of  the  local  Stanton  and  local 
length  Nusselt  numbers.  Comparison  of  the  Stanton  numbers  shows  that  the 
major  effeet  of  shifting  from  the  1-DFM  to  the  2-DFM  is  to  decrease  the 
ealeulated  values  of  the  Stanton  numbers  and  simultaneous1  y to  increase  the 
corresponding  values  of  the  length  Reynolds  number.  In  Fig.  2,  for  a diameter 
Reynolds  number  of  33,000,  the  first  six  values  of  the  Stanton  number  lie  with- 
in an  average  deviation  of  less  than  10  per  cent  from  the  theoretical  solution 
for  tube  flow,  the  remaining  eight  values  fluctuate  about  a mean  value,  and 
the  fluctuations  begin  at  a length  Reynolds  number  based  on  the  2-DFM  of 
600,000.  Similar  fluctuations  are  evident  irt  the  values  of  the  Stanton  number 
based  on  the  1-DFM  and  are  discussed  in  (1).  A similar  process  of  levelling- 
off  of  the  values  of  the  Stanton  number  with  fluctuations  abo’rt  a mean  value  is 
evident  in  Figs.  3 and  4 at  larger  values  of  the  diameter  Reynolds  number. 

Figs.  2 to  7 indieate  that  as  the  value  of  the  diameter  Reynolds  number 
increases  from  33,000  to  115,000  the  values  of  the  Stanton  number  based  on 
the  2-DFM  move  from  below  the  lower  line  for  the  tube  solution  to  almost 
perfeet  agreement  with  the  theoretical  plate-flow  solution,  shown  in  Fig.  7. 

The  discrepancies  between  the  values  of  the  Stanton  number  and  the  tube  solu- 
tion should  not  be  considered  significant  at  present,  sinee  this  tube  solution 
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is  based  on  the  assumption  of  constant  viscosity  and  thermal  conductivity;  work 
is  in  progress  to  obtain  tube  solutions  for  variable  values  of  these  properties. 
The  discrepancies  between  the  values  of  the  Stanton  number  and  the  flat-plate 
solution,  especially  at  the  smaller  values  of  the  diameter  Reynolds  number 
should  not  be  considered  significant  for  two  reasons.  First,  the  plate -solution 
is  based  on  zero  pressure  gradient  but  the  data  were  taken  with  an  adverse 
pressure  gradient  in  the  tube;  note  that  the  best  agreement  between  the  tube- 
flow  data  and  the  flat-plate  solution  occurs  for  the  smallest  value  of  the  adverse 
pressure  gradient  in  the  tube.  Second,  at  the  smallest  values  of  the  diameter 
Reynolds  number,  there  is  some  evidence  of  separation  of  the  laminar  boundary 
layer  so  that  one  would  expect  some  deviations  from  the  theoretical  solution 
for  plate  flow  in  which  separation  effects  were  neglected. 


Nusselt  Number 

The  values  of  the  local  length  Nusselt  number  based  on  the  2-DFM  and 
on  the  1 - DFM  are  also  shown,  in  Figs.  2 co  7.  The  main  reason  for  including 
the  Nusselt  number  is  to  compare  the  value  of  plotting  the  calculations  in 
different  ways.  All  the  conclusions  reached  so  far  with  the  aid  of  the  curves 
of  Stanton  number  are  also  evident  from  examination  of  curves  for  the  Nusselt 
number,  except  that  the  levelling -off  of  the  Stanton  numbers  in  Figs.  2,  3,  and 
4 is  almost  completely  masked  in  the  sharply  rising  values  of  the  Nusselt 
number  at  the  same  length  Reynolds  number.  This  masking  is  a result  of  using 
the  tube  length  as  a multiplicative  factor  in  the  Nusselt  number. 

Transition 

The  transition  from  a laminar  to  a turbulent  boundary  layer  is  clearly 
evident  in  Figs.  5,  6,  and  7,  occurring  as  a very  sharp  rise  in  value  of  the 
Stanton  or  Nusselt  number.  On  the  basis  of  the  2-DFM,  this  sharp  rise  occurs 
at  values  of  the  length  Reynolds  number  of  4,000,000,  4,500,000  and  4,500,000 
in  Figs.  5,  6,  and  7,  respectively.  These  values  of  the  length  Reynolds  number 
at  the  start  of  transition  are  in  good  accord  with  values  given  in  the  literature 
for  transition  of  a laminar  boundary  layer  on  a flat  plate  with  zero  pressure 
gradient. 


Results  for  Test  Combination  C - Effect  of  Shorter  Test  Section 

The  results  of  using  the  2-DFM  for  calculation  of  one  run  with  test  com- 
bination C are  shown  in  Fig.  8.  The  value  of  the  inlet  diameter  Reynolds 
number  is  127,000,  which  corresponds  to  the  highest  vaiue  for  test  combination 
D shown  in  Fig.  7.  Hence  the  effect  of  a shorter  length  of  test  section  may  be 
obtained  by  comparison  of  these  two  figures.  Fig.  8 shows  that  the  value  of 
the  Mach  number,  based  on  the  2-DFM,  remains  nearly  constant.  The  boundary 
layer  thickness  ratio  increases  from  10  per  cent  at  inlet  to  about  20  per  cent 
at  exit,  thus  duplicating  the  behaviour  of  the  boundary-layer  thickness  ratio  of 
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test  combination  D in  Fig.  7.  This  ratio  in  Fig.  8 also  agrees  well  with  th..»t 
computed  by  the  Howarth  formula. 

The  values  of  the  local  Stanton  and  local  Nusselt  numbers  based  on  the 
2-DFM,  shown  in  Fig.  8,  agree  within  a few  per  cent  with  the  values  based  on 
the  theoretical  predictions  for  plate  flow,  in  the  region  where  a laminar 
boundary  layer  is  present.  The  value  of  the  length  Reynolds  number  at 
transition  from  a laminar  to  a turbulent  boundary  layer  lies  between  3,000,000 
and  4,000,000;  this  agrees  well  with  the  values  observed  for  transition  in  the 
longer  test  section  of  test  combination  D.  The  main  effect  of  the  shorter 
length  of  test  combination  C is  that  the  rise  in  values  of  the  Stanton  or 
Nusselt  numbers  during  transition  is  never  as  spectacular  as  that  found  in 
test  combination  D. 


PREDICTION  OF  DIABATIC  SUPERSONIC  FLOW  IN  A TUBE 

The  2-DFM  provides  a means  of  predicting  the  flow  behaviour  for  dia- 
batic  supersonic  flow  in  the  entrance  region  of  a tube  if  a laminar  boundary 
layer  exists.  In  essence,  thi  3 u'c  .hod  o f prediction  was  confirmed  when  test 
combination  D was  designed  to  be  50  diameters  long  instead  of  30  for  the 
previous  test  combination  C.  The  analysis  for  prediction  for  diabatic  flow’ 
is  obtained  by  rearrangement  of  the  basic  equations  given  above  but  will  be 
omitted  here  for  lack  of  space.  The  analysis  baseci  on  the  2-DFM  is  also  being 
used  to  design  new  test  sections  for  adiabatic  and  diabatic  flow  with  an  inlet 
Mach  number  near  5. 


CONCLUSIONS 

An  approximate  and  simplified  two-dimensional  flow  model  for  the  en- 
trance flow  region  of  a tube  yields  results  which  are  in  agreement  with  theo- 
retical predictions  for  tube  flow  and  for  plate  flow  with  zero  pressur  e gradient. 
This  flow  model  yields  a better  understanding  of  the  phenomena  which  occur 
in  diabatic  supersonic  flow  of  air  in  a round  tube  than  that  obtained  with  the 
simple  one -dimensional  flow  model  and  eliminates  the  inadequacies  of  the  latter 
model. 

1 he  Mach  numbers  based  on  the  2-DFM  are  significantly  larger  and  much 
more  nearly  constant  along  the  test  section  than  those  based  on  the  1-DFM. 

The  thickness  of  the  laminar  boundary  layer  based  on  the  2-DFM  in- 
creases siowly  from  tube  inlet  to  exit  but  does  not  fill  the  tube  cross  section 
at  exit.  These  values  of  the  thickness  are  in  good  agreement  with  those  calcu- 
lated from  Howarth' s results. 

The  values  of  the  local  Stanton  number  based  on  the  2-DFM  agree  best 
with  the  theoretical  predictions  for  tube  flow  and  deviate  most  from  those  for 
plate  flow  at  the  lowest  value  of  the  diameter  Reynolds  number,  namely,  33,000. 
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But  the  reverse  is  true  at  the  highest  value  of  the  diameter  Reynolds  number, 
namely,  127,000  for  test  combination  C,  and  115,000  for  test  combination  D. 
The  agreement  between  rneas'.r  ed  tube  values  and  predictions  for  plate  flow 
for  zero  pressure  gradient  for  acre  last  two  values  of  the  diameter  Reynolds 
number  is  within  a few  per  cent.  The  values  of  the  local  Nusselt  number 
based  on  the  2-DFM  behave  similar  to  the  Stanton  numbers  in  comparisons 
with  the  theoretical  predictions  for  tube  flow  and  plate  flow  but  are  of  less 
value  in  understanding  and  interpreting  the  phenomena. 

The  effect  of  varying  the  length  of  test  section  from  30  to  50  diameters 
on  the  results  based  on  the  2-DFM  is  negligible,  except  that  the  shorter  length 
does  not  delineate  transition  as  clearly. 

The  transition  from  a laminar  to  a turbulent  boundary  layer  for  diabatic 
supersonic  flow  in  a tube  is  clearly  evident  in  the  sudden  and  very  sharp  rise 
of  the  vshife  of  the  local  Stanton  number.  The  values  of  the  length  Reynolds 
number  at  the  start  of  the  transition  are  in  excellent  agreement  with  similar 
data  for  a laminar  boundary  layer  on  a flat  plate. 
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APPENDIX 

CALCULATED  RESULTS  BASED  ON  TWO-DIMENSIONAL  FLOW  MODEL 


This  appendix  contains  the  results  based  on  the  2-DFM  for  seven  runs, 
of  which  six  were  made  with  test  combination  D and  one  with  C.  The  summary 
of  the  large  amount  of  computation  based  on  the  2-DFM  is  given  in  Table  1. 

Th^  values  given  in  Table  1 are  based  on  calculations  made  with  five  or  six 
significant  figures  throughout.  A sample  calculation  is  given  here,  in  part, 
co  indicate  the  method,  for  Run  No.  B-l.  The  original  data  for  Run  No.  B-l 
are  taken  from  Table  5 of  reference  (l). 


The  values  are  summarized  as  follows: 


C/ 


w 

* 


w 


w 


46.73  lbm/sec  ft 


(Re*)  = 92,765 

D s 


--  0.967 

= 0.01837  lbm/sec 
= 0.01777  lbm/SwC 


A =■  0.1978  in.“  = 0.00)373  ft' 

From  Tabic  30  of  tlie  "Gas  Tables"  (7),  corresponding  to  the  value  of 

(p/p  .)  of  0.04601  for  the  fourth  station,  it  is  found  that 
oi 


M 


= ?.  655 


l /T  . = 0.4149 

c oi 


= (T  / T .)  T . = 0.4149  (570.88)  = 236.9°  F abs 
c ox  oi 


For  the  core, 

P 

c 

V 


p (p/p  .)  14.479  x 0.04001  x 144 

oi  oi 


RT  53.342  x236.9 

c 


0. 0076 45  lbm/ ft 


3 


M 


y gkJRT  ' = 2.6  55  J 32. 174  x 1.400  x 53.342  x 236.9  = 2003  ft/sec 


= p V = 0.007645  x 2003  = 15.31  lbm/(sec  ft") 

c c ‘ 


For  the  boundary  layer, 
V, 


= V 1 2 = 1002  ft /sec 
c 


From  Equations  (1)  (2),  (3),  (7),  (12),  (13),  (14),  and  (15) 


n=j  - 1 2 

a .12  + S s q + wc  (T  . - V / 2gc  ) - (pV  c A)/R 
\1  tEB*  n p oi  b p h 

A = 

c 


n=l 


c (T  . - vA/2gc  ) - (pV  c )/R 


c o oi 


b p 


q . /2  + > ' q = q /2  + q = C. 00 1482  + 0. 01 1486  = 0.012968  Btu/sec 

j “i  n 4 ■’ 


(T  . - V”  / 2gc  ) = (570.88  - 1002  x 1002/2  x 32.174  x 778.3  x 0.239S) 
oi  b p 

= 487.4  °F  abs 


wc  (T  . - V~/2gc  ) = 0.01777  x 0.2399  x 487.4  = 2. 077  Btu/scc 
p oi  b p 

G c (T  . - V2/ 2gc  ) = '5.31  x 0.2399  x 487.4  * 1791  Btu/sec  ft2 
c p oi  n p 

r 

pV,  c /R  = 14.479  x 0.0460!  x 144  x 1002  x.  0.2389/53.342  = 435.1  Btu/sec  ft ‘ 
b p 

nV  c A/R  = 435.1  x 0.001373  = 0.5976  Btu/sec 
o p 

. 0.012968  + 2.077  - 0.5976  , . , ,2 

Ac  ■ 1791--  435.! 001101  '* 

Ab  = A - Ac  = 0.001373  - 0.001101  = 0.000272  ft2 

w = G A = 15.31  x 0.001101  = 0.01687  lbm/sec 

c c c 

w = w - w = 0.01777  - 0.01687  = 0.000902  lbm/sec 
b c 

5/a  = 1 - (1/a)  y7  Ajx  = 1 - (24/0.5018)  y 0.001 1 01  / tt  = 0.1045 
From  Equations  (3)  and  (7) 


T,  = pV,  A / Rw.  = 14.579x0.04601  x 144x  1002x0.000272/53.342x0.00090 
b b b b 


= 547.3  °F  abs 
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p . (p/p  .)  14.579  x 0.4G01  x 144  _ 

C1  Ol  ) 

= 0.003308  lbm/ft 


RT. 


53.342  x 547.3 


G,  - = 0.0C3308  x 1C02  = 3.314  lbm/sec  ft 

b b b 

The  remaining  quantities  are  calculated  as  follows: 

T = (T  IT  . ) T . - 0.940  x 570.88  = 536.6  °F  abs 
aw  aw  oi  oi 

h - q/A'  (t  - t ) = 974.7/(672.2  - 536.6)  = 7.19  Btu/hr  °F  ft" 
w aw 

hD/A  = 7.19  x 0.5018/12  x 3600  x 0.1942  x 10“5  = 43.0 
c 

hL/A  5 (hD/A  ) (L/D)  - 43.0  x 4.484  - 193 
c c 


h/c  G 5 7. 19/0. 2393  x 15  1 x 3600  = 5.44  x 10 

P c 


-4 


G D (L/D)  15.31  x 0.5018  x 4.484 

He  5 — = — = 463,000 

C ^c  12  x 0.620  x 1C“5 

The  calculations  are  continued  by  repeating  the  above  stops  for  Station 
numbers  5,  6,  7,  etc. 
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C!GURE  6 


REVISED  BISTRIBUTION  LIST  FOR  UNCLASSIFIED 
TECHNICAL  REPOETS  ISSUED  UNDER 
CONTRACT  N5ori-C7805  PROJECT  NR  061-028 


Chid*  of  Naval  Research 
Department  of  t-ne  Navy 
Washington  25,  D*  C, 
Attm  Code  U38 


(2 


' / 


Commanding  Officer 
Naval  Ordnance  Laboratory 
White  Oak,  Silver  Spring,  Md. 
Attn*  Dr,  R,  K.  Lobb 


(1) 


Commanding  Officer 
Office  of  Naval  Research 
Branch  Office 
150  Causeway  Street 
Boston  10,  Massachusetts 

Command! nr  Officer 
Office  of  Naval  Research 
Branch  Office 
31*6  Broadway- 
New  York  13,  New  Tork 


(2) 


v-u 


Commanding  Officer 
Office  of  Naval  Research 
Branch  Office 
Tenth  Floor 

The  John  Crerar  Library  Bldg, 

86  East  Randolph  Stroet 
Chicago  1,  Illinois  (l) 

Conmanding  Officer 

Office  of  Naval  Research 

Branch  Office 

1000  Geary  Street 

San  Francisco  i^  California  (l) 


Commanding  Officer 
Office  of  Naval  Research 
Branch  Office 
1030  E.  Green  Street 
Pasadena  1,  California 


(1) 


Commanding  Officer 
Office  of  Naval  Research 
Navy  #100,  FPO 
New  York,  New  York 

Chiefs  Bureau  of  Aeronaut!  as 
Department  of  the  Navy 
Washington  25,  D*  C, 

Attn:  Research  Division 

Chief,  Bureau  of  Ordnance 
Department  of  the  Navy 
Washington  25,  D,  0, 

Attn?  Cede  Re  9a 


[2) 


(1) 


(1) 


Chief,  Bureau  of  Ships 
Department  of  the  Navy 
Washington  25,  D.  C. 

Attn*  Research  Divirion 

Command lag  Officer  and  Director 
David  Taylor-  Model  Basin 
Washington  7,  D,  C. 

Attn*  Aeromechanics  Division 
Hydrodynamics  Division 


(1) 


(1) 

(1) 


Directorate  of  Intelligence 
Headquarters,  U.S.  Air  Force 
National  Defense  Building 
Washington  25,  D.  C, 

Attn*  Doc,  and  Dissem.  Br.(AFOIR-DD)  (l) 

Armed  Services  Technical 
Information  Agency 
Document  Service  Center 
Knott  Building 

Dayton  2,  Ohio  (5) 

Commanding  General,  Hdqs, 

Air  Research  and  Development  Command 

Office  of  Scientific  Research 

?.  0.  Box  1395 

Baltimore,  18,  Maryland 

Attn*  Fluid  Mechanics  Div.  (l) 

Director,  Research  and  Development 
Department  of  the  A,-^ 

National  Defense  Building 
Washington  25,  D«  C,  (l) 

Commanding  General 
Office  of  Ordnance  Research 
Department  of  the  Army 
Washington  25,  D,  C,  ( 


1} 


Ballistic  Research  Laboratories 
Department  of  the  Army 
Aberdeen  Proving  Ground,  Maryland 
Attr.s  Hr,  J,  Sternberg  (l) 
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Director  of  Rasearch 
National  Advisory  Committee 
for  Aeronautic* 

172U  F Street,  Northwest 
Washington  25,  D.  C.  (l) 

Director 

National  Advisory  Committee 
for  Aeronautic* 

Langley  Aeronautical  Laboratory 
Langley  Field,  Virginia  (l) 

Director 

National  Advisory  Committee 
for  Aeronautic* 

Langley  Aeronautical  Laborato:  y 

Langley  Field,  Virginia 

Attn:  L'r.  A.  Bueemann  (l) 

IB  rector 

National  Advisory  Committee 
for  Aeronautics 
Ames  Aeronautical  Laboratory 
Moffett  Field,  California  (l) 

Director 

National  Advisory  Committee 

for  Aeronautics 

Ames  Aeronautical  Laboratory 

Moffett  Field,  California 

A.ttn:  Dr.  D.  R.  Chapman  (1) 

Director 

National  .Advisory  Committee 
for  /-ronautics 

Lewis  Flight  Propulsion  Laboratory 

21000  Brookpark  Road 

Cleveland  11,  Ohio  (l) 

JlT-wcuor 

National  Advisory  Committee 
for  Aeronautic* 

Lewie  Flight  Propulsion  Laboratory 

21000  Brookpark  Road 

Cleveland  11,  Ohio 

Attn*  Dr*  J.  C.  Eward  (l) 

The  Johns  Hopkins  University 
Department  of  Aeronautical  Engineering 
Baltimore  18,  Maryland 
Attn:  Dr.  F.  H.  Clauser  (l) 

The  Johns  Hopkins  University 
Aoplied  Physics  Laboratory 
8621  Georgia  Avenue 

Ml?  u) 


University  of  Michigan 

Department  of  Aeronautical  Engineering 


East  Engineering  Building 
Ann  Arbor,  Michigan 

Attn:  Drs  Arnold  Kuathe  (l) 

Prof.  E.  W.  Conlon  (1) 

California  Institute  of  Technology 
Guggenheim  Aeronautical  Laboratory 
Pasadena  U,  California 
-Attn:  Dr.  Hans  W.  Liepmaan  (l) 

Prof.  Lester  Lee*  (l) 

Dr.  P,  A.  L&gerstrom  (1) 


Director 

National  Bureau  of  Standards 
Department  of  Commerce 
Washington  25,  D.  C. 

Attn:  Dr.  G.  B.  Schubauer  (l) 

Polytechnic  Institute  of  Brooklyn 
Department  of  Aeronautical  Engineering 
and  Applied  Mechanics 
99  Livingston  Street 
Brooklyn  2,  New  Tork 

Attn:  Dr.  A.  Ferri  (]) 

Brown  University 

Division  of  Engineering 

Providence  12,  Rhode  Island 

Attn?  Dr.  ?.  F.  Maeder  (l) 

California  Institute  of  Technology 

Jet.  Propulsion  Laboratory 

Pasadena  California 

Attn:  Mr.  F,  E,  Goddard  (l) 

University  of  California 

Engineering  Research  Projects 

Berkeley  h,  California 

Ar.t.n:  Dr.  S.  A,  Schaaf  (1) 

University  of  California 

Department  of  Engineering 

Los  Angeles  2L»>  California 

Attn:  Dean  L.M.K.  Boelter  (1) 

Cornell  Aeronautical  Laboratory 
]iU55  Gsnessee  Street 
Buffalo  21,  New  York 

Attn;  Mr.  A.  H.  Flax  (l) 

Cornell  University 

Graduate  School  of  Aeronautical 

Engineering 

Ithaca,  Nev?  York 

Attn:  Dr.  W.  R.  Sears  (i; 
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University  of  Delaware 
ueparwionv  ox  wn«iiaC«i. 

Newark,  Delaware  . . 

Attn:  Dr,  Kurt  Wohl 

Harvard  University 

Department  of  Applied  Physics 

and  Engineering  Science 

Cambridge  38,  Massachusetts 

Attn  Prof-  H#  W.  Emmons  P/ 

Massachusetts  Institute  of  Technology 
Guided  Missiles  Couanittee 
Cambridge  39,  Massachusetts 
Attn:  Supervisor,  Tech.  Reports  vlJ 

Massachusetts  Institute  of  Technology 
Department  of  Aeronautical  Engineering 
Cambridge  39,  Massachusetts 
Attn:  Dr.  M.  Finston  ^ 

Massachusetts  Institute  of  Technology 
Department  of  Mechanical  Engineering 
Cambridge  39.,  Massachusetts 
Attn i Prof*  J,  H.  Keenan  (1J 

Prof.  A.  H,  Shapiro  (l) 

Princeton  University 

Dapartrocirt  of*  Aaro nautical  Engin® citing 

Princeton,  New  Jersey 

Attn:  Prof.  S,  M.  Bogdonoff  (1^ 


Purdna  University 

School  of  Mechanical  Engineering 

Lafayette,  Indiana 

Attn:  Dr,  N.  J.  Zucrov  (l) 

University  of  Texas 
Defense  Research  Laboratory 
300  East  2Uth  Str&et 

Austin,  Texas 

Attn:  Dr,  M,  J.  Thompson  (l) 

Rensselaer  Polytechnic  Institute 
Department  of  Aeronautical  Engineering 
Troy,  New  York 

Attn:  Dr,  R.  P.  Harrington  (l) 

University  ctV  Maryland 

Institute  for  Fluid  Dynamics 

and  Applied  Mathematics 

College  Parx,  Maryland 

Attn:  Dr,  S,  I,  Pai  (l) 

Case  Institute  of  Technology 
Department  of  Mechanical  Engineering 
Cleveland,  Ohio 

Attn:  Prof.  G.  Kuerti  (l) 

Stanford  University 

Department  of  Mechanical  Engineering 

Stanfo  rd,  Calif  omis 

Attn:  Prof.  A.  L,  London  (l) 
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Because  of  our  limited  supply,  you  are  requested  to  return  this  copy  WHEN  IT  HAS  SERVED 
YOUR  PURPOSE  so  that  it  may  be  made  available  to  other  requesters.  Your  cooperation 
will  be  appreciated. 
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NOTICE:  WHEN  GOVERNMENT  OR  OTHER  DRAWINGS.  SPECIFICATIONS  OR  OTHER  DATA 
ARE  USED  FOR  ANY  PURPOSE  OTHER  THAN  IN  CONNECTION  WITH  A DEFINITELY  RELATED 
GOVERNMENT  PROCUREMENT  OPERATION,  THE  U.  3.  GOVERNMENT  THEREBY  INCURS 
NO  RESPONSIBILITY,  NOR  ANY  OBLIGATION  WHATSOEVER:  AND  THE  FACT  THAT  THE 
GOVERNMENT  MAY  HAVE  FORMULATED,  FURNISHED.  OR  EN  ANY  WAY  SUPPLIED  THE 
SAID  DRAWINGS,  SPECIFICATIONS,  OR  OTHER  DATA  IS  NOT  TO  BE  REGARDED  BY 
IMPLICATION  OR  OTHERWISE  AS  IN  ANY  MANNER  LICENSING  THE  HOLDER  OR  ANY  OTHER 
PERSON  OR  CORPORATION,  OR  CONVEYING  ANY  RIGHTS  OR  PERMISSION  TO  MANUFACTURE, 
USE  OR  SELL  ANY  PATENTED  INVENTION  THAT  MAY  IN  ANY  WAY  BE  RELATED  THERETO. 
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